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Abstract 

Quasipatterns (two-dimensional patterns that are quasiperiodic in any 
spatial direction) remain one of the outstanding problems of pattern for- 
mation. As with problems involving quasiperiodicity, there is a small 
divisor problem. In this paper, we consider 8-fold, 10-fold, 12-fold, and 
higher order quasipattern solutions of the Swift-Hohenberg equation. We 
prove that a formal solution, given by a divergent series, may be used to 
build a smooth quasiperiodic function which is an approximate solution 
of the pattern-forming PDE up to an exponentially small error. 

Keywords: bifurcations, quasipattern, small divisors, Gevrey series 
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1 Introduction 

Quasipatterns remain one of the outstanding problems of pattern formation. 
These are two-dimensional patterns that have no translation symmetries and 
are quasiperiodic in any spatial direction (see figure [1]) . In spite of the lack 
of translation symmetry (in contrast to periodic patterns), the spatial Fourier 
transforms of quasipatterns have discrete rotational order (most often, 8, 10 
or 12-fold). Quasipatterns were first discovered in nonlinear pattern-forming 
systems in the Faraday wave experiment [10, 14], in which a layer of fluid is 
subjected to vertical oscillations. Since their discovery, they have also been 
found in nonlinear optical systems [19], shaken convection [29,33] and in liquid 
crystals [26] , as well as being investigated in detail in large aspect ratio Faraday 
wave experiments [1,4,5,25]. 
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Figure 1: Example 8- fold quasipattern. This is an approximate solution of 
the Swift-Hohenberg equation ([T]) with /i = 0.1, computed by using Newton 
iteration to find an equilibrium solution of the PDE truncated to wavenumbers 
satisfying |k| < \/5 and to the quasilattice r27. 

In many of these experiments, the domain is large compared to the size of 
the pattern, and the boundaries appear to have little effect. Furthermore, the 
pattern is usually formed in two directions (x and y), while the third direction 
(z) plays little role. Mathematical models of the experiments are therefore often 
posed with two unbounded directions, and the basic symmetry of the problem is 
E{2), the Euclidean group of rotations, translations and reflections of the {x,y) 
plane. 

The mathematical basis for understanding the formation of periodic patterns 
is well founded in equivariant bifurcation theory [16]. With spatially periodic 
patterns, the pattern-forming problem (usually a PDE) is posed in a periodic 
spatial domain instead of the infinite plane. Spatially periodic patterns have 
Fourier expansions with wavevectors that live on a lattice. There is a parame- 
ter /i in the PDE, and at the point of onset of the pattern-forming instability 
(/X = 0), the primary modes have zero growth rate and all other modes on the 
lattice have negative growth rates that are bounded away from zero. In this case, 
the infinite-dimensional PDE can be reduced rigorously to a finite-dimensional 
set of equations for the amplitudes of the primary modes [8, 9, 17, 21, 32], and 
existence of periodic patterns as solutions of the pattern-forming PDE can be 
proved. The coefhcients of leading order terms in these amplitude equations can 
be calculated and the values of these coefficients determine how the amplitude 
of the pattern depends on the parameter fi, and which of the regular patterns 
that fit into the periodic domain are stable. Due to symmetries, the solutions of 
the PDE are in general expressed as power series in y^/I, which can be computed, 
and which have a non-zero radius of convergence. 

In contrast, quasipatterns do not fit into any spatially periodic domain and 
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have Fourier expansions with wavevectors that hve on a quasilattice (defined 
below). At the onset of pattern formation, the primary modes have zero growth 
rate but there are other modes on the quasilattice that have growth rates arbi- 
trarily close to zero, and techniques that are used for periodic patterns cannot 
be applied. These small growth rates appear as small divisors, as seen below, 
and correspond at criticality (/x = 0) to the fact that for the linearized opera- 
tor at the origin (denoted —Cq below), the eigenvalue is not isolated in the 
spectrum. 

If weakly nonlinear theory is applied in this case without regard to its va- 
lidity, this results in a divergent power series [30], and this approach does not 
lead to a convincing argument for the existence of quasipattern solutions of the 
pattern-forming problem. 

This paper is primarily concerned with proving the existence of quasipatterns 
as approximate steady solutions of the simplest pattern-forming PDE, the Swift- 
Hohenberg equation: 

^^^,U-{l + AfU-U' (1) 

where U{x,y,t) is real and fi is a parameter. We do not prove the existence 
of quasipatterns as exact steady solutions of the PDE. We are not concerned 
with the stability of these quasipatterns: in fact, they are almost certainly un- 
stable in the Swift -Hohenberg equation. Stability of a pattern depends on the 
coefficients in the amplitude equations (as computed using weakly nonlinear 
theory). In the Faraday wave experiment, and in more general parametrically 
forced pattern forming problems, resonant mode interactions have been identi- 
fied as the primary mechanism for the stabilisation of quasipatterns and other 
complex patterns (see [31] and references therein). These mode interactions are 
not present in the Swift-Hohenberg equation, though their presence would not 
significantly alter our results. 

In many situations involving a combination of nonlinearity and quasiperiod- 
icity, small divisors can be handled using hard implicit function theorems [13], 
of which the KAM theorem is an example. Unfortunately, there is as yet no 
successful existence proof for quasipatterns using this approach, although these 
ideas have been applied successfully to a range of small-divisor problems arising 
in other types of PDEs [12,22,23]. There are also alternative approaches to 
describing quasicrystals based on Penrose tilings and on projections of high- 
dimensional regular lattices onto low-dimensional spaces [24] . 

We take a different approach in this paper: we show how the divergent power 
series that is generated by the naive application of weakly nonlinear theory can 
be used to generate a smooth quasiperiodic function that (a) shares the same 
asymptotic expansion as the naive divergent series, and (b) satisfies the PDE ([1]) 
with an exponentially small error as ^ tends to 0. This approach is based on 
summation techniques for divergent power series: see [2,7,28] for other examples. 
In order to make the paper self-contained, we put in Appendices some proofs of 
useful results, even though they are "known" . 



3 



In section[21 we define the quasilattice and derive Diophantine bounds for the 
small divisors that will arise in the nonlinear problem, for Q-fold quasilattices: 
Lemma [2.11 extends the results of [30] covering the cases Q = 8, 10, 12, to any 
even Q > 8. We then compute in section [3] (following [30]) the power series 
in ^JJi for a formal Q-quasipattern solution U of the Swift-Hohenberg equation, 
where /i is the bifurcation parameter in the PDE. 

In section 131 we define an appropriate function space Tia'- each term in the 
formal power series U is in this space. In section [5l we prove (Theorem 15. ip 
bounds on the norm of each term in the formal power series solution of the PDE. 
In the Q-fold case, the norm of the /z"+2 term in the power series for the 
quasipattern is bounded by a constant times ii'"(n!)'*', where if is a constant 
and Z + 1 is the order of the algebraic number lo = 2 cos(27r/(5), which is also 
half of Euler's Totient function (p(Q) (Z = 1 for Q = 8, 10 and 12, / = 2 for 
Q = 14 and 18, I = 3 for Q = 16, 20, 24, 30, . . . ). This result was announced 
in [30] for Q < 12, and is extended here to Q > 14. With a bound that grows 
in this way with n, the power series is Gevrey-Al, taking values in a space of 
Q-fold quasiperiodic functions. 

In section [6l for convenience, we consider the cases Q = 8, 10 and 12. We 
introduce a small parameter ^ related to the bifurcation parameter by C = -yZ/I, 
so that the norm of the ^^"+2 term in the power series for U is also bounded 
by a constant times K"{n\)'^ < if"(4n!). We use the Borel transform U of the 
formal solution U: the ^''"+2 term in the power series for U is the ^''"+2 term 
in the power series for U divided by (An + 2)!. With this definition, U is an 
analytic function of ( in the disk |^| < if^^/^, and for each C in this disk, U 
is a Q-fold quasiperiodic function of {x, y) in the space Tis- Of course the new 
function U does not satisfy the original PDE, but we prove that it satisfies a 
transformed PDE (Theorem 16. 2p . 

The next stage would be to invert the Borel transform: however, the usual 
inverse Borel transform is a line integral (related to the Laplace transform) 
taking ( from to oo, and U is only an analytic function of ( for ^ in a disk. 
If the definition of U could be extended to a line in the complex ( plane, the 
inverse Borel transform would provide a quasiperiodic solution of the PDE - 
this remains an open problem. 

Since the full inverse Borel transform cannot be used, in section [71 we use 
a truncated integral to define jJ{i'). This involves integrating C, along a line 
segment inside the disk where U is analytic, weighted by an exponential that 
decays rapidly as v ^ 0. We show that t/(j^) and U{fi) have the same power 
series expansion when we set ly = ^fji, but unlike J7, lJ(y) is a C°° function 
of in a neighbourhood of 0, taking values in Jis- In other words, is a 

Q-fold quasiperiodic function of (x, y) for small enough /i. This function is not 
an exact solution of the Swift-Hohenberg PDE, but we show in Theorem 17.21 
that the residual, when fj{ij}^^) is substituted into the PDE, is exponentially 
small as /i ^ 0. Finally, in the last section [8l we show that by taking as initial 
data the above approximate solution, the time dependent solution U{t) stays 
exponentially close to the approximate solution for a long time, of the order 
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O(l///1+1/40. 

In conclusion, we have shown that, for any even Q > 8, the divergent power 
series (7 (/it) generated by the naive apphcation of weakly nonlinear theory can 
be used to find a smooth Q-fold quasiperiodic function U{^^/^^) that shares the 
same asymptotic expansion as J7, and that satisfies the PDE with an exponen- 
tially small error. 

This technique does not prove the existence of a quasiperiodic solution of 
the PDE. However, this is a first step towards an existence proof for quasiperi- 
odic solutions of PDEs like H]). In particular, we may hope to use C7 as a 
starting point for the Newton iteration process that would form part of an ex- 
istence proof using the Nash-Moser theorem. As an aside, ordinary numerical 
Newton iteration succeeds in finding an approximate solution of the truncated 
PDE for values of /i where the formal power series has already diverged, as in 
figure [3l 

An analogous result may be proved for example in the Rayleigh-Benard con- 
vection problem (see [20]), using the fact that the dispersion equation possesses 
the same property as in the present model: at the critical value of the parame- 
ter there is a circle of critical wavevectors in the plane. The method might also 
extend to the case of the Faraday wave experiment by considering fixed points 
of a stroboscopic map. 

In the present work we consider quasilattices generated by regularly spaced 
wavevectors on the unit circle, and solutions invariant under 2t: /Q rotations. It 
might be worth studying the case of solutions having less symmetry on the same 
quasilattice, or quasilattices (still dense in the plane) generated by wavevectors 
that are irregularly spaced. 

Acknowledgments: We are grateful for useful discussions with Sylvie Ben- 
zoni, W. Crawley-Boevey, Andre Galligo, Ian Melbourne, Jonathan Partington, 
David Sauzin and Gene Wayne. We are also grateful to the Isaac Newton In- 
stitute for Mathematical Sciences, where some of this work was carried out. 



2 Small divisors: Quasilattices and Diophantine 
bounds 

Let Q e N (Q > 8) be the order of a quasipattern and define wavevectors 

kj = (cos (2^^) ' sin (^2^^^^^ , j = 1, 2, . . . , Q 

(see figure [5^). We define the quasilattice F C to be the set of points spanned 
by integer combinations of the form 

Q 

km = ^ Tijkj, where m = (toi, m2, . . . , mg) g N*^. (2) 
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(a) 



(b) 



(c) 






Figure 2: Example quasilattice with Q = 8, after [30]. (a) The 8 wavevectors 
with |k| = 1 that form the basis of the quasilattice. (b,c) The truncated quasi- 
lattices Fg and T27- The small dots mark the positions of combinations of up 
to 9 or 27 of the 8 basis vectors on the unit circle. Note how the density of 
points increases with A''k. 



The set F is dense in R^. 

We are interested in real functions U{x.) that are linear combinations of 
Fourier modes e^^'^, with x e and k e F. If [/(x) is to be a real function, we 
need Q to be even, with kj and — kj in F, hence the quasilattice F is symmetric 
with respect to the origin. 

In the calculations that follow, we will require Diophantine bounds on the 
magnitude of the small divisors. We see below that the small divisors are 
||k|^ — l|, for k e F. To compute the required lower bound, we start with 



^ 2mj,mj^ cos(ji - ^2)^0 + ^ 



TO" 



l<jl<j2<Q 

where ^0 = Stt/Q. Let us define 



i<j<Q 



u) = 2 cos ^0 



We now show how |k„ 



l| can be expressed as a polynomial in co. 



First, we can express 2 cosp^o as a polynomial in w, for 1 < p < Q — 1: 



2 cosi 



■ pui- 



P-2 



P(p-3) ^,p-4 



with integer coefficients which only depend on Q (easy proof by induction), and 
the leading coefficient being 1, and since cos(p + Q/2)9q = —cosOq, this leads 
to 

|k™|'= «X> ?;eZ,r = 0,l,...,g/2-l, (3) 

0<r<Q/2-l 

where the integers q'^. are quadratic forms of m. 

Next, we use the property that u is an algebraic integer, since it is the sum 
of two algebraic integers e*^° + e^*^"^)*^". More precisely, a; is a root of the 
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(minimal) polynomial P{uj) with integer coefficients, with leading coefficient 
equal to 1, and which is of degree (p{Q)/2 := I + 1, where (p{Q) is Euler's 
Totient function [3], the number of positive integers j < Q such that j and Q 
are relatively prime. For example, <p(14) = 6 since the 6 numbers 1, 3, 5, 9, 11 
and 13 have no factors in common with 14, and so Z + 1 = 3 in the case Q = 14. 
In the cases Q = 8, 10 and 12, the irrational numbers uj = 2 cosOo are V2, ^^./^ 
and \/3: these are quadratic algebraic numbers + 1 = 2), while for Q = 14, uj 
is cubic. 

Finally, dividing ([3|) by P{lo) we obtain a remainder of degree / such that 

\km\^-l=qa+u}qi + ---+Jqi (4) 

where qo + 1 and qj, j = 1, . . . , Z are integer- valued quadratic forms of m. 

Define |m| = ruj, then, for a given wavevector k e F, we define the order 
TVk of k by 

iVk-min{|m|;k = k„„k„eF}. (5) 

The reason for this is that, for a given k, there is an infinite set of m's satisfying 
k = km- For example, we could increase mj and mj+Q/2 by 1: this increases |m| 
by 2 but does not change km- Whenever solutions are computed numerically, 
it is necessary to use only a finite number of Fourier modes, so we define the 
truncated quasilattice F^y to be: 

Fjv - {k e F : A^k < iV} . (6) 

Figure mb,c) shows the truncated quasilattices Fg and F27 in the case Q = 8. 
For example, we have in the case Q = 8: 

4 

|km|^ — ^'^rn'j~ + V2{m[ni2 + rn2r)i'^ + rn'^rn'^ — rn'^m'i) , (7) 

4 

iVk = (8) 

where m'^ — rrij — ^j+Q ji ■ More generally we have 

Q/2 

i=i 

The above inequality can occur strictly (for example) in the case Q = 12, 
because only 4 of the 12 vectors kj are rationally independent in this case. 
More generally only tp{Q) vectors kj are rationally independent [34]. 

Now, the quantity in (|4]), \qo +ujqi + ■ ■ ■ + oj''qi\, may be as small as we want 
for good choices of large integers qj , and we need to have a lower bound when 
this is different from 0. 
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In [30], it was proved that in the cases Q = 8, 10 and 12, there is a constant 
c > such that 

||kp - l| > for any k e r with |k| 7^ 1. (9) 

The proof relies on the fact that for quadratic algebraic numbers, there exists 
C > such that 

C 

\p -m\> — 

holds for any (p, g) G Z^, g 7^ [18]. Now using the fact that q is quadratic in 
m (see ([7])) we have 

q < QNl (10) 

from which ([9]) can be deduced. 

The Diophantine bound ^ may be extended to any even Q > 8, and there 
exists c > depending only on Q, such that for any k e F, with |k| 7^ 1, we 
have 



c 

To show this, we use the following known result (see [11]) proved in Appendix El 



1 > (11) 



Lemma 2.1 Let uj he an algebraic number of order I + 1, that is, a solution 
of P{lij) — where P is a polynomial of degree I + 1 with integer coefficients, 
that is irreducible on Q. Then, there exists a constant C > such that for any 
Q ~ (it), qi, ■ ■ ■ , qi) G Z'^"'^\{0}, the following estimate 

C 

\qo + qiUJ + q2Uj'^ + ■ ■ ■ + qiuJ^\ > Yjl (12) 

holds, where |q| = Eo<j<i l^jl- 

In the general case, by choosing m such that |m| ~ iVk„, the estimate (|10p 
is replaced by 

|q| < c{Q)Nl 

where c{Q) depends only on Q. Then estimate ([TT]) is satisfied by taking 

- ^ 

It remains to show that |k|^ 7^ 1 for all k e F, apart from k ~ ki, . . . , kg. 
This is solved by denoting ^ = e'^" , the Qth primitive root of unity, and relat- 
ing kj_|_i to C"' , and k^ to ™i+iC"' - We then use the Kronecker- Weber 
theorem which says that "every abelian extension of Q is cyclotomic" [34] . This 
implies that if "T-jC"* (which is an algebraic integer) has modulus 1, then it is 
necessarily a root of unity. Knowing that the dimension of the Q- vector space 
spanned by the C"' is v{Q), this implies that this root of unity is one of the 
j = l,...,Q. 
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3 Formal power series computation 



Let us consider the steady Swift-Hohenberg equation 

{l + AfU-^iU + U^ = (13) 

where we look for a Q-fold quasiperiodic function [/ of x e K^, defined by 
Fourier coefficients C/k on a quasilattice F as defined above. Let us rewrite (|13p 
in the form 

CoU = nU -U^ (14) 

where 

We write formaUy 



C/(x) = ^ke' 



ker 

the meaning of this sum being given in section 31 We seek a solution of 
bifurcating from the origin when = 0, that is invariant under rotations by 
27r/Q. First we look for a formal solution in the form of a power series of an 
amplitude. More precisely we look for the series 



C/(x,m) 



./fEA*"t^^"H-) (15) 



as a formal solution of (fT3)) . where all factors t/*^"-' are invariant under rotations 
by 2tt/Q of the plane. The coefficient /3 will be given by fixing 
At order 0{y/\li\) in (fT3|) we have 



CqU'-"^ = (16) 

and we choose the solution 

Q 

t/(o) ^^g»k, x^ (17) 

which is invariant under rotations by 2^ /Q and defined up to a factor which we 
take equal to 1. 

In writing in this way, we have made use of the fact that the only 
solutions k G F of |k| = 1 are k = ki, . . . , kg (see discussion at the end of sec- 
tion[2l). This implies that the kernel of is only one-dimensional if restricted to 
functions invariant under rotations by 2-k/Q, this kernel being spanned by U''^\ 

At order 0(1^^/^) we have 

/:oC/(i^ =C/(") -/3-\C/("))^ (18) 

We need to impose a solvability condition, namely that the coefficients of e^^^'^, 
for j = 1, . . . , Q on the right hand side of this equation must be zero. Because of 
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the invariance under rotations by 2tt/Q, it is sufficient to cancel the coefficient 

/3 = 3(Q-1)>0, (19) 

and U^^^ is known up to an element /3^^''[/*^°) in ker£oi which is determined at 
the next step: 

[/(I) = C7(i)+/3(i) C7(i)= "ke^*""' (20) 

ker,Ark=3 

3 

ask, = -1/64, a2k,+k, = ~ (i _ |2k^- + k;P)2 ' kj + k/ 7^ 0, 

= -(i-ik,+'L,+k.p)2' 

kj + k, ^ 0, kj + k, ^0, k, + k, ^ 0, 

where f/^^-' has no component on e^^^'^. 
Order in ^ leads for n > 2 to 

= /J- 1 ^ C/Wc/(') (21) 

fc,i,r>0 

For n — 2, the solvability condition on the right hand side gives /3(^\ and f/(2) 
is then determined up to /j'-^-'C/^'^^ . Indeed we obtain on the right hand side 

{/(I) _ ^t7(l)[/(0)2 ^ ^(1) ^(l)[/(0) _ i^(l)[/(0)3 _ 3 ^(l)j^(0)2 

= -2/3(i)c/(") + C/(i) - ^L^(i)[/(o)2 _ ^£ot/''\ (22) 

P P 

where we used the fact that the component of [/(o)-^ on e^^^'^ is /3 (see (IT51) ). 
Hence 2/3^^ is the coefhcient of e'^^ '' in -S/J^^t/^i^C/^''^^^ and since all coeffi- 
cients of U^^^ are negative, we find /J^^-* > 0. We obtain in the same way the 
coefficients /J^^-^^C/'"' of [/("-i) in using the solvability condition on the right 
hand side of (1211. 



Small divisor problem. It is clear that we can continue to compute this 
expansion as far as we wish, where at each step we use the formal inverse of Cq 
on the complement of the kernel. However, applying Cq^ to eJ^'^ introduces a 
factor 

1 

(l-|kP)^' 

which may be very large for combinations k = km with large m, since points 
km of the quasilattice T sit as close as we want to the unit circle. This is a small 
divisor problem and computations indicate that the series (jlSp seems to diverge 
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0.5 




Figure 3: Amplitude A^^^ of the quasipattern, as a function of ^ and of N, with 
Q = 8, iV = 1, 3, 9 and 27, and scaled so that A^^^ = y/JI. Increasing the order of 
the truncation leads to divergence for smaller values of fi. The squares represent 
amplitudes computed by solving the PDE by Newton iteration, truncated to the 
quasilattice T27 (Ni^ < 27) and restricted to wavevectors with |k| < VS- Note 
that for fj, = 0.1, the Newton iteration succeeds in finding an equilibrium solution 
of the PDE, while the formal power series has diverged. The spatial form of the 
solution with /i = 0.1 is shown in figure [T] 

numerically [30] . We illustrate this in figure [31 plotting the amplitude A^^^ 
against ^, where 




and the norm || • \ \s and the projection operator Pq are defined below: A^^^ is 
essentially the magnitude of the coefficient of e^^^ '^ as a function of /i and of N, 
the maximum order of wavevectors included in the truncated power series. 

However, we prove in section [5] that in all cases we can control the diver- 
gence of the terms of the series (|15p . and obtain a Gevrey estimate ||?7''"-'||s < 
7ii'"(n!)'", where the norm || • ||s is defined below. 

Remark 3.1 For Q — A or 6, there is no small divisor problem since T is a 
periodic lattice, and the only points in T that lie in a small neighborhood of the 
unit circle are — 1, • . • , Q}- 
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4 Function spaces 

We characterise the functions of interest by their Fourier coefficients on the 
quasilattice T generated by the Q unit vectors k^: 

C/(x) = ^ke^" " 
ker 

Recah that for each k e F, there exists a vector m e N'^ such that k = km = 
J2f=i ^j^i ^"^d '^^ choose m such that |m| = iVk as defined in ([S]). We 
have the following properties, proved in Appendix [Bl 

Lemma 4.1 (i) We have the following inequalities: 

A^k+k' < A^k + A^k' , A_k - A^k, (23) 

|k| < Ak. (24) 

(ii) We have the following estimate of the numbers of vectors k having a 
given N\^: 

card{k : A^k = A} < ci(Q)A^'5/2"i (25) 
where Ci{Q) only depends on Q. 

Define now the space of functions 

ns^\u = J2Uke^''-- ||C/||? = 5](l + Ak2)^|C/kp<ool, (26) 
I ker ker J 

which becomes a Hilbert space with the scalar product 

{W,V),^Y.^l + NkyWkV^. (27) 
ker 

In the sequel we need the following lemma, proved in Appendix [Cl 

Lemma 4.2 The space Tig is a Banach algebra for s > (3/4. In particular there 
exists Cs > such that 

||C/l-||.<c.||C/|M|F||,. (28) 
For £ > and s > I + (5/4, Tis is continuously embedded into . 

From now on, all inner products are s unless otherwise stated, so that we 
can remove the s subscripts throughout in scalar products. 

We will also use the orthogonal projection on ker£o: for any U ^Tis, let 

and we denote by Qq the orthogonal projection: 

Qo =I-Po, 

which consists in suppressing the Fourier components of e'*'^ '', j — 1, ... ,(5- 
The norm of the linear operator Qo is 1 in all spaces TLs- 
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5 Gevrey estimates 

In this section we prove rigorously a Gevrey estimate of in (fT5|) . The 

estimate for Q = 8, 10 and 12 {I = 1) was announced in [30]. Recall that a 
formal power series X^^o^"^" Gevrey-k [15], where fc is a positive integer, 
if there are constants S > and i^T > such that 

\un\<SK"-{nl)^ Vn>0. (29) 

Theorem 5.1 For any even Q > 8, assume that s > (3/4. Then there exist 
positive numbers K{Q, c, s) and 5{Q, s) such that there exists a unique formal 
solution U (fi) of U3\) . under the form of a power series in /i^^^, all factors [7*^") 
belonging to Tis, and which satisfies 




U(n) ^ +[/("), (t/("\e*^-"), =0, j = l,...,Q, 

l^(«)| < jii:«(„!)4i^ 

where I = ■^^{Q) — 1 is the integer defined in Lemma \2.1\ From the above 
inequalities, it follows that 

||f/(")||, <7i^"(n!)4', „>0, 

where 7 is related to S, Q and s only. 

Remark 5.2 The above Theorem claims that the series U in powers of is 
Gevrey-2l taking its values in Tig. 

Remark 5.3 In the cases when Q — A or 6, the pattern is periodic, and the 
above series may be built in the same way, leading to a series which is convergent 
for 11 < fiQ where /ip > 0. This results simply, via the Lyapunov-Schmidt 
method, from the implicit function theorem in its analytic version. The values 
of fiQ for Q — 2, 4 or 6 are estimated in [30]. 

Proof. We choose s > Q/A since Lemma |4?2] insures that Tis is then a Banach 
algebra. We notice that 

||e»k,.x||^ ^2^/2^ 

and 

= 2^/2yQ. (31) 

We also have [3^°^ = 1 and [/(°) = 0. Now we notice from ^ that for |k| 7^ 1 
we have 

9 TV*' 

||k|^-l| (32) 
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which controls the unboundedness of the pseudo- inverse £q ^ (inverse of £o 
restricted to the orthogonal complement of its kernel). Indeed Cq^ is bounded 
from Hs to Tis-Ai- 

Remark 5.4 We may notice that the set of eigenvalues of Cq is dense in the 
positive real line, which constitutes the spectrum. Hence is not isolated in the 
spectrum of £o- This explains why the pseudo-inverse of Cq on the complement 
of its kernel, is unbounded and satisfies (see iS^) ): 

||A)"'QoC/||.-4i < ^\\U\U, for any U G Tls. 

The basic observation here is that the factor [/'^"^ that multiplies /i" has a 
finite Fourier expansion in e^^'^, with k = Y^f^i "n^j^j, < 2n + 1, hence 

Ny^ < 2n+l. Since for C/^^-' we have |m| = 3 in all k^'s, equation (|18p leads to 

We set ^ ^ 

lj(n) ^ ^(")c/(0) + [/(")^ ^{n) ^ Qg[/(")^ (34) 

and replacing this decomposition in (|2ip we obtain, by taking the scalar product 
with e'^^^'^ 

^("-l)2-_i/3C/(«-l)f/(0)2^g.ki.x\ i / ^ik)JJ^l)JJ(r)^^^W,■A 

0<k,l,r<n-2 

where we have used {U^"\e^^^'^) — We^^^'^Wj. = 2". Next, we use 

(3[/(«-l)[/(o)2^g^kl.x^ = (3^7(°)^e*l•'') + (3L^("-i)L/(°)2^e*i-^) 

= 3/3/3(""^)2^ + (3C/("-i)C/(")2^ e^*'!-^), 

and we are led to solve with respect to L/'^"-' the following system for 

n > 2 



rof/H = [/("-I) _ /3-1q„ U^'''^U^^^U^^'\ (35) 

k+l+r=n~l, 
k,l,r>0 

f^(n-i) ^ _L /3^("-i)[/(")2 + ^ C/Wc/Wc/M,e'^''-A.(36) 

^ \ k+l+r=n-l, I 

0<k,l,r<n-2 

Now we make the following recurrence assumption: there exist positive constants 
7]^, S and depending on Q and s, such that 

m^P^ls < liKP{p\Y\ p = 0,l,...,n-l, (37) 
< SKP[p\f\ p=l,...,n-2. 
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These estimates hold for J7^"^ = and for U^^'' provided that and K satisfy 

Putting these together resuhs in 

\\U^P^\l = \\p^P'>U^°^+U^P^\l < (2^/^(5v/Q + 7i) KP{pl f, 

or 

<7XP(p!)4', with 7 = 2^/2^v^ + 7i. (39) 

The resolution of ^ and ^ provides and f/^"^ starting with n = 2. 

A useful lemma is the following, proved in Appendix ID] 

Lemma 5.5 The following estimates hold true for I > 1: 
n3,„ = (k^-l^-rT < 4(n!)^', n > 1 



k-\-l-\-r—n 



3,n 



^ (fc!Z!r!)4' < 10((n- 1)!)*', 7i > 2. 

0<fc,i,r<n-l 

Thanks to Lemma [5751 and the estimate for in ([55)1 . we observe that 



r) 



fc+/+r— n— 1 
0</cJ,r<n-2 



From this it follows that 



< 10c273x"-i((n-2)!)4'. 



I/?^"-'' I < ^T^^""'((" - 1)0'' {37i2^Q + 107^} , 
and the recurrence assumption is realized if 

9 

{37i2"Q + 1073} < 5 



3(Q- 1)21+^/2 
holds. Now we have, still by using Lemma [5751 



(40) 



(41) 



< i^"(n!)4'(2 + -)4'-i- 



71 + y 7^ 



The factor (2n + l)'*'/c2 here comes from the pseudo-inverse of Cq acting on 
functions containing modes of order up to 2?i + 1 . The recurrence assumption 
is realized if 



7i + -jl^ \ < liK. 



(42) 
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We now must choose 7]^, S and K in such a way as to satisfy the three 
conditions ((38)) . (|40)) and (jlS]) . Indeed, we may choose such that 

and replacing this value in ([39]) and (|40p . then ([40]) is satisfied as soon as 

-3 



2 3(Q- 1)2-/2-1 / Q-i 



holds. Then, choosing K such that 



,34'/ 2«+ic2g\ 1 34'-i22'^c4Q5/2 
A = max < 1 ' » 1 



c2 V 5{Q-1)J ' 6 c2(g- 1)2 



allows to satisfy ([38]) and (|42|). 

We conclude that the bounds on ||C/(")|U and |^("^| in Theorem [Q hold, 
and that ([39|l . which holds for < p < n — 1, also holds for p = n, and so 



||C/(")||, < 7i^"(n!)4', n>l. 
This ends the proof of Theorem 15.11 ■ 

6 Borel transform of the formal solution 

In this and subsequent sections, we consider the cases with ^ = 1 (Q = 8, 10 
and 12) and set 

Remark 6.1 In the general case, we should set = fi^^^K 

The formal expansion ([50]) becomes, after incorporating (3^^^^ into U^'^\ 

U ^C^Y.C*''u'-"\ (43) 

n>0 

and we have the estimate 

\\U^"^\\s < jK"{n\)^ < 7ii:"(47i!). 

Thus the formal power series is a Gevrey-1 series in (. 

Let us now consider the new function C i— > U{C,), taking its values in Tig, 
defined by 



4n+2 



^ (An + 2)1 



n>0 



Indeed, by construction, this function is analytic in the disc |C| < ^ = K ^1^^ 
with values in the Hilbert space TYs and invariant under rotations of angle 2t: /Q. 
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The mapping U ^ U, where we divide the coefficient of C" by n\, is the Borel 
transform [6] apphed to the series U . Since U satisfies a Gevrey-1 estimate, the 
Borel transform U is analytic in a disc. 

We now need to show that this function C/(C) is solution of a certain par- 
tial differential equation. Let us recall a simple property of Gevrey-1 series. 
Consider two scalar Gevrey-1 series u and v 

n>l n>l 
\Un\ < CiK^nl, \Vn\ < C2K^nl, 



then we have 



{uv)n = E UkVn-k, 

l<k<n-l 

\{uv)n\ < CiC2K{'n\, 



as this results from Appendix [Dl by using the following inequality for n > 3 
^ E fc!(n-fc)!<l + 2(i + ... + -l^)<n, 

^ ' l<fe<n-l 

which shows that in our case we can multiply two Gevrey-1 series with coef- 
ficients belonging to Tig (the factor C1C2 is then multiplied by Cg) and obtain 
a new Gevrey-1 series with coefficients in Hs. It is then classical that we can 
write ^ 

*gU *gU (44) 
where the convolution product^ written as *Gi is well defined by 

UkVn-k 



n>l l<fe<ri-l 

and satisfies 

{u *G v) — (uv). 

This convolution product is easily extended for two functions /(C) and g{C), 
analytic in the disc |C| < K^^, and with no zero order term, by 

(/*5)(C) = E E A5n-.^^^^C". (45) 

ri>ll<fc<n-l 

It is clear that for / = u, and g — v we have 

/ * 5 = (w *c {)) = (uv). 
Since we have (l44l). it is clear from (|2T]) that we have 



((l + A)2C/)(x,C) = (l + A)2[/(x,C). 
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Now let us define a bounded linear operator K, as follows: for any function 
C ^ ^(C) analytic in the disc |C| < iff^, taking values in Tis, canceling for 
C = 0, and satisfying 

n>l 

we define 

n>l ^ '' 

It is then clear for V ~U that 

n>0 ^ 

and we see that 

dl{K.U) ^ U. 

We now claim the following: 

Theorem 6.2 The Borel transform J7(x, C) o/ the Gevrey solution found in 
Theorem \5.1\ for I = 1 is the unique solution, analytic in the disc \(^\ < K"^^^, 
cancelling for ( = 0, and taking values in Tig invariant under rotations of angle 
27r/(5, of the equation 

(1 + l^fV -K.V + V *V *V ^Q. (46) 

Proof. We assume ^ = 1 in what follows. The changes needed for larger Vs are 
left to the reader. Let us look for a solution V in the form 

n>l 

where Vn G TLg is invariant under rotations of angle 2'it/Q. Then defining a 
formal series 

C/ = E C"^"' = n!K, 

n>l 

it is clear that U satisfies formally 

(1 + A)2;7-C^f7 + C/^ = 0, 
and by identifying powers of C^: 

CoUi = 0, 
C0U2 = 0, 
JO-oUs + U! = 0, 
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which leads to ?7i = because of the last equation where the solvability condi- 
tion cannot be satisfied. Then we have 

[/i=0, £ot/, =0,j = 2,3,4,5, 

and 

CoUq ~U2 + = 0. 

We observe that U2 and Uq satisfy the equations verified by /3~^^^[/^°^ and 
^-i/2j|(i) ^ggg (jYgj) )^ This is indeed the only solution invariant under rotations 
oi 2Tr/Q. Hence 

Now at order ('^ we get 

CoUv - f/g + WiUs = 
and since C/3 = CU^'^\ where C is a constant, the solvability condition gives 

hence C = and C/3 = 0. It is the same for U4 — U5 = 0, and we obtain 
^oUr = CqUs = CoUg 0. Then the computation of higher orders is exactly 
as the one for the computation of C/*-"', since the cubic term cancels if the sum 
of the 3 indices p in Up is not 2 mod 4. Coming back to the definition of 
Un = nlVn, it is then clear that Theorem l6.2l is proved. ■ 

7 Truncated Laplace transform 

Let us take K' > Ki and define a linear mapping C/ 1— > C/ in the set of Gevrey-1 
series taking values in Tis 

U{v)^- e-iu{C)dC, (47) 
^ Jo 

where C/(C) is the Borel transform of U as defined above, which is analytic in 
the disc |C| < 1/Ki. The function i/ C7(i/) is a truncated Laplace transform 
of the Borel transform of U. 

Remark 7.1 // C/(C) could be shown to be analytic on a line in the complex 
^ plane extending to 00, instead of just in a disk, then the Laplace transform 
in |^7| j would be the inverse Borel transform, and would provide a quasiperiodic 
solution of U3\) in Hs- 



19 



It is clear that U{u) is a C°° function of in a neighborhood of 0, taking its 
values in Hs, as this results from 



U{v) = / e-^U{iyz)dz 
Jo 

and from the dominated convergence theorem. Moreover (/{u) and U{iJ,) have 
the same asymptotic expansion in powers on v, when we set n = v^^^, as this 
results from 

' e ~^dC^ i/"-e K'^ [ — • • • 



1 K'll K'^-^ln- 1)! 

(48) 

It is also clear that in a little disc near the origin 

U = U, 

but this does not imply that U = U since U is not a function, being defined 
as a formal series of z/^, and an asymptotic expansion does not define a unique 
function. The real question is whether or not U is solution of in Tis- 
By construction, we know that the Gevrey-1 expansion of 

y(^l/4) ^. (1 + A)2[7(^l/4) _ ;,[7(^l/4) + [7(^1/4)3 

in powers of //^Z"* is identically 0, but we don't know whether this function 
(smooth in /i^^"*), which is in 7is-4, is indeed 0. In fact we have the following 

Theorem 7.2 For any even Q > 8, take s > Q/i. Then, I = {l/2)tp{Q) - 1 

being defined by Lemma \2.1{ the quasiperiodic function [/(/i^/"*') £ Tig, with 

s > Q/4, defined from the series found in Theorem \5.1[ is solution of the Swift- 

1 

Hohenberg PDE \13\) up to an exponentially small term bounded by C{K')e k'^^/"-' 
in Hs-i, for any K' > K^/^K 

Proof. The result of the Theorem follows directly from two elementary lemmas 
IE. II and IE. 21 on Gevrey-1 series shown in Appendix [El and which may be un- 
derstood in the function space Tis instead of C. Indeed, for Z = 1 this gives an 
estimate of the difference beween F (/i^/^) and the truncated Laplace transform 
of the left hand side of equation (^5]) (which is then 0), taking into account of 

(1 + ^fu{^i^/') = 4u r ^"^(1 + ^fuiOdC, 



//1/4 ; 
' JO 

which holds in Hs-i- Using Rcmark l6.1| the extension to larger Z's is left to the 
reader. ■ 
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8 On the initial value problem 

Once we know an approximate solution U of the steady PDE (|13p. a natural 
question is: let us start at time t — Q with U\t=a — U, what can we say about the 
solution U{t) of the initial value problem, for t > 0? Let us give the following 
partial answer to this question: 

Lemma 8.1 Assume Q > S and s > (3/4 and consider the solution U{t) of the 
initial value problem |ip with U\t=o — U € Ti.s+4, where U is given by Theorem 
\7.2\ Then there are a and C" > such that the estimate 

\\U{t)~U\\s<C'e'^ 

holds for < t < ^1+1/41 , where c is the same as in Theorem \7.S\ 

Proof. We can replace s in Theorem 17.21 bv s + 4, hence we have 

with C and c > such that 

\\U\\s+4<C^,, \\R\\s<Ce-^. 

Let us introduce the semi-group e^'-'°*,t > 0, defined for any U £ TCs, s > 0, by 

(e-^ott/)k = e-(i-|k|')'*[/k. 

This semi-group is strongly continuous in Tig, and bounded by 1. Now defining 
W{t) = U{t) - U, we have in Hs 

W{t) = [ e-^°^*-^'>{^lW{T)~3U^W{T)-3UW{T)'^ -W{Tf}dT + 
Jo 

+ [ e-^°^'-^'^RdT. (49) 

We know that W{0) = 0, and by standard arguments the solution of the initial 
value problem exists at least on a finite interval [0, T) in Hs- Let us give a more 
precise estimate on W{t) for a part of the interval of time where < 
C\^fji for a certain Ci > 0. A simple estimate on p9)) leads to 

IIW^WII. < r72||W^(r)|Mr + iCe"^, 
Jo 

with 

72-(l + 3C2-|-3CCi + Ci2)^. 
Then solving this inequality by Gronwall, we obtain 

\\Wms< (e-^^*-!) 

72 

which leads directly to the result of the Lemma. ■ 
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A Proof of Lemma 12.11 



We give below an elementary proof of Lenima l2.1l 

The polynomial P being irreducible on Q of degree I + 1 and the polynomial 
Q defined by 

Q{x) = ^ qjx\ 
o<j<i 

being of degree I, then by the Bezout Theorem there exist two polynomials A{x) 
of degree I — 1 and B{x) of degree Z, with coefficients in Q such that 



A{x)P{x) + B{x)Q{x) = 1. 



(50) 



Defining coefficients Pj , < j < Z + 1 , , < j < ^ — 1 and 6^ , < j < Z of 
polynomials P, A and B, the identity ([50]) becomes a linear system of 2Z + 1 
equations, of the form 

MX = ^0, (51) 



where the unknown is X with 



X = 



M = 



/ Pi+i 
Pi 



Pi 
Po 


V 





pl+1 



Po 



( \ 

a.j-2 

bi 
bi-i 

V ^0 J 






Pi+i 
Pi 
Pi-i 

Pl-2 











Po 



qi 

qi-i 
qi-2 

qo 










V 1 J 



qi 

qi-i 





qi 



qo 




\ 





qi 

qi-i 



qo J 



The (2Z + 1) X (21 + 1) matrix M has integer coefficients and is invertible (other- 
wise it would contradict the Bezout Theorem). Hence its determinant is integer 
valued and is an homogeneous polynomial of degree I + 1 in q = (qo, qi) . 
We may invert the system (fSTj) by Cramer's formulas and we observe that the 
coefficients bj are rational numbers, with a common denominator of degree I + 1 
in q and with a numerator of degree I only (we replace in the determinant one 
column containing the g^'s by ^g). It results that the polynomial B[x) is the 
ratio of a polynomial with integer coefficients Bq of degree I in q, with an integer 
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d, homogeneous polynomial of q of degree / + 1 and which is different from 
(detM ^ 0). Now taking x = uj in ^ leads to 



and since d>l and the coefficients of Bq are bounded by C"|q|', this completes 
the proof of Lemma 12.11 

B Proof of Lemma 14.11 

Assertion (ii) follows from the fact that we can group the coefficients rrij — 
"^j+Qji — i ^-nd since in the Q/2~ dimensional space of {m'p j = 1, . . . , Q/2\ 
the set I = is a union of 2*^/^ simplexes of area of order 0(A'5/^^^). 

To prove the part (i) (I23p we observe that 

Q 

A^k+i = min{|m + 7i|; k + 1 = ^(mj + ?ij)kj} 

Q Q 

< min{|m|; k = rrij'kj} + min{|n|; 1 = njkj} 

< Ak + Ai, 

where 

Q Q 
TVk = min ^^m^kj, iVi = min ^^n^kj. 
^'^"^ 3=1 

We notice that Nq = 0, and A_k ~ A^k (each m'j for k is just the opposite for 
— k); we deduce that inequality (^5]) may be strict, since 

= A^o = Ak-k < A^-k + A^k = 2A^k. 
The last inequality ([M)) is easily deduced from 

Q 

k = ^mjkj 

where {nij} gives precisely the "norm" A'k; which implies (since m| — 1) 

Q 

|k|<^|m,|=A^k, 
i=i 

and the Lemma is proved. 
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C Proof of Lemma 14.21 

Let u e Tls, then by Cauchy-Schwarz inequality in P{T) (F is countable) we 
have 











ker 





Now by 



ker " ^ 

we have the following estimate 



Q/2-1 



ker 



iieN 



„ik-x 



which is bounded when s > Q/4,. Hence for s > Q/A the series X^ker^ke' 
converges absolutely and represents a continuous quasiperiodic function, the 
norm (uniform norm) of which being bounded as soon as the norm in Tis is 
bounded. We may proceed in the same way for the derivatives in using (|24|) . 
and show that the series 



ker 



is absolutely convergent for s > Q/4+L This ends the proof of the last assertion 
of the Lemma. Let us now prove the first assertion which is necessary for our 
nonlinear problem. 

First step: We first use the following inequality due to (I23p 



(1 + iV^+k')'/' < 2'"' {(1 + ^^k)'^' + (1 + ^k')'^'} 



valid for any s > 1, because of (1231) and a simple convexity argument (this 
inequality is in fact valid for s > 0). Then the following decomposition holds 

2 



E 

K 



E 

k+k'=K 



with 



Si 



E 

K 

- E 



E "kl'k' 
k+k'=K 



2^s 



K 



E 

k+k'=K 



(1 + ^k) 



(1 + ^k')^ 



For symmetry reasons in the space (k, k'), it is then sufficient to estimate Si. 
Let us split the bracket in the sum Si into two terms: a sum S[ containing 
(k, k') such that 

iVk < 3iVk', 
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and a sum S*" containing (k, k') such that Ny^ > 3iVk'. Hence we have now 

^1 < 2{S[ + S'l) 



with 



K 



E 

k+k'=K, 



s'^ = E 



K 



E 

k+k'=K, 
JVk>3Ark' 



-UkWk' 



2ns 



To estimate S[ we use which gives A'k < 4A^k', hence 



1 16 

< 



and, in using again Cauchy-Schwarz 



'K 



|UkWk' 



k+k'=K. 

JVk<3A'k/ 



It results that 



k+k'=K. 



' (l+iV2-)«/2 



< 



(l+iV2.)-/2 



s[<\\ur^j\v\\'n^J2a 

which, for s > Q/A leads to 



K 



S[<C\\u\\j,J\v\\j^^. 

Second step: We now find a bound for S", which is more technical, since we 
split this sum into packets of increasing lengths. 
Let us define 

ApU = ^ uue'^'^, A_iu = Mo- 

2P<Nk<2P+^ 

It is clear that for s > Q/A (the series is absolutely convergent) 

oo 

u = ApU. 

p=-l 
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Moreover, it is clear from the definition that the norm of u £ Hg is equivalent 
to 

/ oo \l/2 

To estimate the sum S'{, we notice that in the product uv the terms ApuAqV 
only take into account the wavevectors k and k' such that 

2P<7Vk<2P+\ 2«<iVk/<2«+\ 7Vk>37Vk'. 
TVk' < 2P, 2«+i < TVk, 



This implies 
hence in S" 



ApuAqV = 0, for p < q. 
Now, we use (for the sum in S'^) 



K 



k+k;'=K, 
JVk>3iVk/ 



2 \s 



and the right hand side is the square of the norm of the product uv computed 
on terms such that A'^k > SiVk', k + k' = K. We now use the equivalent norm 
defined above with the decomposition in packets, hence 



p>0 \g=-l / 



Let us define Sp-iv = J2q=-i ^9^' then we have 



hence by Cauchy-Schwarz 



Now 



Aj Y Sp-ivApi 



< 



\p=j-i 



Y 2'P^\\Aj{Sp.,vApu)\\l 
p=j-i 



\\Sp-ivApu\\l = ^ I Y "kt^k' 

K k+k'=K,0<iVk,<2P<Ark<2P+i 
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and a classical computation (convolution * P) using Cauchy-Schwarz gives 

K k+k'=K K I k+k'=K k' J 

< (EK'i)(EEi^K-kikkp) 



k K 

2 



which leads to 



\Sp^ivApu\\l < \\Apu\\l (E bk' 



and since the series^ |wk'| < c||t;||-H^ for s > (3/4, as shown at the beginning of 
the proof of Lemma 14. 2( we have 



^vApu\\'^<C\\Apu\\l\\v\\^^^. 

Finally, we obtain 

i+i i+i 
E 2'P^\\A,iSp-^vApu)\\l < E 2""' \\Sp^ivApu\\l 
P=j-i p=j-i 

< C'\\v\\'^^ E 2''1l^p"llo> 
P=j-i 

and 

<C"||HI?,, E 2''1l^f"llo- 



Aj Sp-ivApuj 



p=j-i 



hence 



S'l < 3C"\\v\\'^^ E 2'''1|Ap«llo 
< C'lll^ll^Jbllk 

and Lemma 14.21 is proved. 

D Proof of Lemma 15.51 



Let us define the two sums 

n' 



ri-l 

J2iki{n-ky.) 



k=l 
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we have already 

n2,o = 1, n2,i = 2, n2,2 = (2 + ^)(2!)4', 

nb = 1, n^,3 = 2(2!)4', 

which shows that Il2,n < (2 + for n = 0, 1, 2, and / > 1. Now we have 

for n > 2 

n2,„+i n2,„ _ g/ fc!(n-fc)! y' // n + l-fc y ^| ^ 



((n+l)!)4' (n!)4' ri\ J [\ n+1 

2 2 2^' 



(n+l)4' (n(n + l))4'' 

and since n^' — (n + 1)*' + 2^'~^ < for n > 1 the above right hand side terms 
are negative. It results that for n > 2 

/ (n+l)! ^' 

i-l2,n+l S j J-J-n,2, 

hence 

n2,n< (2 + ^)(n!)4', n>0. (52) 

In the same way 

n'2,„+i _ n'2,n ^ / fc!(r7.-fc)! \ r / n+l-fc y _ , \ ^ ^ 

{{n-iyr n\ ) \\ n + l ) 

hence for n > 2 

TT' TT' o4i 

- ((n- 1)!)4' + 



and 



n2,„ ^ o4z/ 1 1 



< 24M- — _ + ... + _ +n' 



((n-l!)4' - V(n-l)*' ^24';^"^'^ 



Finally 



n'a.,, < 3((n - 1!)*' for n > 2. (53) 
Consider now n3^„ defined by 

n3,n= E (fci^w)"- 

fe+i+r— n 
fe,;,r>0 
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We already have 

n3,o = l, n3,i=3, n3,2 = (3+^)(2!r' <4(2!f', 
In splitting the sum we obtain easily for n > 3 



n3,„ = n2,„ + (n!)'" + ^(r!)4'n2,„- 



< (3 + ^)(n!r' + (2 + ^)n',,„ 

< + ^ + 3(2+1)-!) 
- (3 + ^ + |)(^!r'<4(n!)4'- 



Hence 



n3,„ < 4(n!)4' (54) 
holds for any n > 0. Consider now Ilg „ defined for n > 2 by 



We already have 



k+l+r=n 
0<k,l,r<n-l 



n'3,2 = 1, 

and for n > 3, we obtain in the same way 



Kn = n^,„ + E(H)4'n2.„_. 

< (^3 + 3(2 + l))(n-l!)4' 

< 10(n-l!)4'. (55) 
Hence, with estimates ([S3]) and ([55]) . Lemma [S751 is proved. 

E Lemmas on Gevrey-1 series 

Below we give elementary proofs of two useful lemmas. The interested reader 
will find more general results in [28] and [27]. 

In the following we denote by Ck' the linear operator defined for analytic 
functions v on the disc {jzj < l/^i} by 



1 



l/A" 



{Ck'v){v) = -1 e--v{z)dz, K' > K 
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We also use the notations 



||w||o,i^' = sup \v{z)\, \\v\\iji, = sup \v'{z)\, 
zeio,i/K') zeio,i/K') 

and when w(0) — 0, we notice that (integrating by parts for the second estimate) 

\iCK'v){iy)\ < \\v\\o.K', (56) 



{CK'v){,y) 



l/K' 



e "v'{z)dz 



< e ff'^llul 



0,K' 



Then we have the following Lemmas giving estimates of the commutator of 
Ck' o S (where B is the Borel transform) with the multiplication by ly^ and with 
the mapping u u'^ in the space of Gevrey series. 

Lemma E.l Assume that u{iy) is a Gevrey-1 series, with uq = 0, then for 
V < l/K' 



For any given Gevrey-1 series u, with uq — 0, there is C{K') > such that for 

V < vq{K') we have the estimate 

(Ck'W^) {y) - {CK'uf {v)\ < C{K')e-ih, K' > K^. 

Lemma E.2 Assume that u^iy) is a Gevrey-l series, with uq — 0, then for 

V < l/K' there exists C{K') such that 

\{CK'JCu){y) - v^Ck'u\ < CiK')\\u\\o,K'e-i^. 
Proof of Lemma lE.li From the identity 



^ z™ ^{z — Zi — Z2)' 

{k~ iy.{m- i)m 



dz2 dzi 



(fc + m + 0!' 



from the definition (|15|) of the convolution product, and from the analyticity of 
u in the disc {|2:| < 1/Ki}, we have 







u' {zx)u' {z2)u{z — Zi — Z2)dz2 I dzi ) dz. 



/o 

By Fubini's theorem and a simple change of variables, we obtain 



u' {zi)u' {z2)u{z'i)dzidz2dz^ (57) 
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where Dk' — {zi, ^2,23 > 0; zi + Z2 + 23 < l/K'}. Now, we have 

e 

and from (|56p we obtain 



'u{zi)u{z2)u{z^)dzidz2dz2,, 



[CK'uf [v) - - 



e " ^ u' {zi)u' {z2)u{z3)dzidz2dz3 



< 



<e '<'''\\u\\l^K'i\\u\\o,K' + i^WuWi.K')- (58) 
Now, we observe that (0, 1/ K'Y\Dk' is such that zi + Z2 + 2:3 > ^/K\ hence 



(0,1/K'Y\D^, 



e ^ " u' {zi)u' {z2)u[z:i)dzidz2dz3 



(59) 

CoUecting (|571) . ([55)1 and ([55)1 the first result of Lemma IE. II is proved. Notice 

that by choosing K" > K' , then for v small enough < e~ . Since 

is chosen arbitrarily larger than Ki, we can assert that u being given, there is 
C{K') such that 

(^Ck,^^ (i/) - {CK'uf H\ < C{K')e-i^, K'>Ki. 



Proof of Lemma IE. 21 By integrating by parts, we obtain 

{£K'ICu){iy) = -e-T^ [{ICu) + v{K.u)' + v^{JCu)" + v^{K.u)"'] \x/k- + 
+v\CK'u){y). 

Hence 

\{Ck'1Cu){v) - v\Ck'u){v)\ < e-i^\\u\\o,K' |^ + ^ + ^ + 
which proves Lemma lE.21 ■ 
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